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ABSTRACT 

We show that a particular many-matrix model gives rise, upon hamiltonian 
reduction, to a multidimensional version of the Calogero-Sutherland model and its 
spin generalizations. Some simple solutions of these models are demonstrated by 
solving the corresponding matrix equations. A connection of this model to the 
dimensional reduction of Yang-Mills theories to + 1-dimensions is pointed out. In 
particular, it is shown that the low-energy dynamics of DO-branes in sectors with 
nontrivial fermion content is that of spin-Calogero particles. 
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The quest for integrable (non-relativistic) particle systems in more than one 
spatial dimension is often frustrating. In general, no nontrivial such systems exist 
(with the exception of some isolated few-body cases), that is, systems with a non- 
quadratic potential and which are not a repackaging of one- dimensional degrees of 
freedom. In one dimension, on the other hand, a whole class of integrable many- 
body systems is known, namely the Calogero model and its various generalizations 
(also known as Calogero-Sutherland-Moser systems) [1-3]. The question then is to 
what extent these models remain solvable, if at all, in higher dimensions. 

A particularly fruitful approach to Calogero-like systems is though matrix mod- 
els, in which the particle positions are regained as the eigenvalues of some appro- 
priate matrix [4,5]. The integrability, as well as the solutions of the equations 
of motion, are simpler to obtain this way. It would seem, then, that this is the 
most promising route to systems of higher dimension. The purpose of this note 
is to show that, indeed, appropriate matrix models give rise, under hamiltonian 
reduction, to multidimensional many-body systems of the Calogero type. These 
matrix models are not in general integrable, and therefore we do not expect the 
corresponding particle systems to be integrable either. The hope is, nevertheless, 
that the (inherently simpler) matrix models dynamics will allow for a better study 
of the dynamics of the particle systems. 

The starting point will be a many-matrix model consisting of d time-dependent 
hermitian N x N matrices Mi, i = 1, . . . d, which we will also represent as a vector 

— # 

of matrices M. The action will be the usual kinetic term for each matrix plus some 
potential invariant under simultaneous unitary conjugation of the matrices. The 
eigenvalues of Mj then will be interpreted as the 2-component of the position vectors 
of iV particles moving in flat d- dimensional space. For this interpretation to be 
possible, however, the matrices must be simultaneously diagonalizable, else there is 
no invariant association of the n-th eigenvalue of the matrices as the coordinates of 
the same particle. The eigenvalues of all Mj can be simultaneously permuted with 
a common unitary transformation, which corresponds to having identical particles. 
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So we write the lagrangian 

} ~ V(M) (1) 



L = tr I \™i + iJ2 l^AMi, Mj] 
{ i ij 



where overdot stands for time derivative. = — Lji is an antisymmetric set of 
d(d—l)/2 hermitian matrices serving as Lagrange multipliers for the commutativity 
constraint between the Mj. The potential V(M) = V^U^ 1 MU) can be any real 
function of the Mj invariant under simultaneous unitary transformations of the 
Mj. The form V = trV(M), where V(x) is some real scalar function on R d , will 
be assumed in what follows, which leads to an external potential V(x) for the 
particles. The harmonic oscillator potential V(x) = ^uj 2 x 2 , corresponding to the 
matrix potential V = ^cj 2 ^jtrM?, is the simplest example. 

For a central potential the model is invariant under SO(d) rotations Rif 



Mi -> RijMj , A i:j -> RikRjiAki (2) 
As a result, there is a conserved angular momentum 

Jij = tr{MiMj - MjMi} (3) 
Time translation invariance implies the conservation of energy 



£ = tr j^M2 + \/({M}) 



(4) 



The invariance of (1) under simultaneous conjugation of all matrices by a time- 
independent unitary matrix implies the existence of a conserved matrix "angular 
momentum" 

K = tJ2[Mi,Mi] (5) 

i 

K is a traceless hermitian matrix. Choosing it to have the "minimal" form where 
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all its eigenvalues are equal but one, that is, 

K mn = i(5 mn - u m u n ) , ^2\u m \ 2 = N (6) 

m 

where u is a fixed iV-vector, will lead to Calogero dynamics for the eigenvalues of 
Mj, just as in the one-dimensional case, 

The equations of motion are 

M i + V i (M)+iY t [A ij ,M j ]=0 (7) 

3 

plus the constraint 

[Mi,M,-]=0 (8) 

The constraint implies that Mj can all be diagonalized with a common time- 
dependent unitary rotation U(t): 

Mi = U~ l XiU , Xi = diag(xi t i, . . . x ijN ) (9) 

In terms of (9) the equations of motion acquire the form 

X i + 2[X i ,A] + [X i ,A]+[[X i ,A],A]+V4X)+i^2[A ij ,X j ]=0 (10) 

3 

where A = UAU^ 1 and A = UU^ 1 is the "gauge potential" generated by the time 
variance of U. 

We now recall that the commutator of a diagonal matrix with any matrix has 
zero diagonal elements, since 

[D, B] mn = (d m - d n )Bij , if D = diag(d\, ...d N ) (11) 

Therefore, isolating the diagonal terms in (10), only the first, fourth and fifth term 
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contribute and we have 



%i,m + ^ ^ 2(xj jjn Xi !n )A mn A nm + V^(:r m ) — (12) 

n 

Plugging the form (6) and (9) in (5), on the other hand we have 

i ^ ^ { x i,m x i,n) A rnn = £(5 mn U m U n ) (13) 



i 



where u = Uu. For m ^ n and m = n we obtain the relations for A mn and u m , 
respectively 

_ _lhj^Ur } _ 2 _-| /-, 4 N 

Plugging these in (12) and calling (xi jTn , . . . Xd, m ) — %m we finally obtain 



- 2f l Xm _ ^4 + = (15) 



which is the equation of motion for the positions of particles x m in an external 
potential V(x) and interacting through a rf-dimensional two-body inverse square 
potential £ 2 /x 2 , that is, a ci-dimensional generalization of the Calogero model. 
The key elements in the derivation are that the Lagrange multiplier term does not 
influence the eigenvalue equations of motion and that the angular part reproduces 
the rotationally invariant d- dimensional inverse-square potential. Note that the 
energy (4) and angular momentum (3) become the corresponding quantities of the 
particle system in the constraint subspace, that is 

E = E l*m + E 75 £ 1 7] 2 + E v ^ ( 16 ) 

Jij = XiXj XjXi (17) 

If the "angular momentum" K is not in the "minimal" form (6), it will enter 
the equations for the eigenvalues in a nontrivial way and will give rise to multidi- 
mensional generalizations of the 'spin-Calogero' model [6-9]. To see this, we point 
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out that the restriction of the Hamiltonian in the constraint subspace [Mi, Mj] = 
takes the form 

H = H \& + E w^r? + E ntn) (is) 

m m^n v ' m 

where K = UKU^ 1 . As usual, K mn Poisson-commute to the SU(N) algebra and 
can be recast into internal degrees of freedom ("spin") for the particles [10]: 

p 

Kmn — Sm^H (19) 
o=l 

To study the matrix equations of motion we specify to the minimum nontrivial 
dimensions d = 2 and to the rotationally invariant harmonic external potential 
V(x) = ^uj 2 x 2 . Defining the non-hermitian matrix M = Mi + 1M2, the equations 
of motion and constraint become 

M + [A, M] + u 2 M = , [M, Mt] = (20) 

while the "angular momentum" K takes the form 

K = i[M\M] (21) 

in the constraint subspace. Solving the two-dimensional Calogero model amounts 
to finding solutions of the above matrix equations for M with the form (6) for K. 

The simplest possible class of solutions is the one with A = 0. It can be 
shown, however, that these solutions correspond to linear motion of the particles 
and the system reduces to the one-dimensional Calogero model. The solution to 
the equations of motion is 

M = Ae iuJt + B^e~ iuJt (22) 
where the matrices A, B, to satisfy the commutativity and "angular momentum" 
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constraints, must obey 



[A,B] = 0, [A,At] = [5, fit] = £ (23) 

with K as in (6). In terms of the new matrices Q = A + P+ and P = iuo(A — P/f) 
(representing the position and velocity matrices at t — 0) relations (23) become 

[Q, Qt] = [P, pt] = , [Q,pt] = -i/f (24) 

This tells us that P and Q, although non-hermitian, can each be diagonalized with 
a unitary rotation, with complex eigenvalues (representing the initial positions 
and velocities of the particles on the complex plane). Choosing a basis where Q 
is diagonal with eigenvalues q m , we deduce from the last relation in (24) that the 
matrix elements of P are 

it 

P-mn = Wm$mn H ; ~(1 $mn) (25) 
Qm 9n 

where we used \u m \ = 1 in the Q-diagonal basis and further chose the phases of 
the states such that u m = 1. From [P,pt] = o now we obtain 

— — — = real , V ^— = real (26) 

Qm ~ Qn {Qm ~ Qk){Q m ~ Qn) 

By using the invariance of the equations by a shift of Q and P by a multiple of the 
unit matrix (which is related to the fact that the center-of-mass motion decouples 
from the relative motion) we can always choose q\ — p\ — 0. Then the first relation 
above implies that all q m are collinear (i.e., g m /g n =real) unless p m = aq m for some 
real a. The second relation, however, is satisfied only if the q m are collinear. By 
the first relation, p m will also be collinear with q m . Therefore, we see that the 
one-dimensional Calogero model is included in the A = sector of the general 
model. 
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The other case in which the equations of motion have an obvious solution is 
when A=constant. Choosing a basis in which A is diagonal, with (real) eigenvalues 
A n , we have 

Mmn + u mn M mn = , u mn = to 2 + A m — \ n (27) 
which has as solutions (we assume tu mn > 0) 

M mn = A mn e^ + Bl m er iUmnt (28) 

The task of finding the most general A mn , B mn which satisfy the commutativ- 
ity and "angular momentum" constraints is not trivial. We demonstrate here a 
particularly simple solution, namely 

Amn = A m 5 mn + CL n 5 m ,n+l i ^ran = B n 5 mn , 

i (29) 
\a n \ — a , co n+ i n+ 2 — Wti, n +i — — 9 

The last constraint for u; n)n +i translates into — 1 algebraic equations for the 
N — 1 variables A n — A n +i (clearly A can be shifted by any multiple of the unit 
matrix). It is obvious that the diagonal part A m , B m represents a general motion 
of the decoupling center of mass. The eigenvalues of the off-diagonal part of A are 
the N-th roots of a\ • • ■ ajv- So the off-diagonal part of M has eigenvalues 

z m = x m + iy m = ae^+^V ? Ur = J2 ^,n+i = oo 12 + (N — l)£/a 2 (30) 

n 

Therefore the relative motion is one in which the particles are regularly positioned 
on a circle of radius \a\ and rotate with constant angular velocity cu r . 

The above model can be generalized to one with unitary matrices. Omitting 
the details of the calculation, we simply state the result. The lagrangian of the 
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model is 



[ i ij ) 

where again V is some real conjugation-invariant potential. The Lagrange multi- 
plier matrices A y - = — Aji are not hermitian, but the constraints arising from the 
variation of A y - and aJ- are compatible (in fact, equivalent) for unitary C/j. The 
eigenvalues of Ui, written as exp(ixi, m / Ri) represent coordinates of particles on a 
d-dimensional torus of radii Ri. Upon choosing the "angular momentum" 

iY J Rl[Ulu i ]=l{l-uv)) (32) 



as before, the x n move like particles on the torus in an external potential V(x n ) 
and interacting through a periodic generalization of the ci-dimensional two-body 
inverse-square potential 



V(x n - x m ) = — — — — (33) 



Similarly, the rf-dimensional generalization of the inverse-sinh model can be ob- 
tained by taking Ri — > iR^. 

It is interesting to note a connection of the models presented here with the 
matrix model obtained as the dimensional reduction of (d+l)-dimensional Yang- 
Mills theory. In the A = gauge, this model is essentially the kinetic part of the 
above model plus a potential of the form 

tr^^Mj-] 2 (34) 

ij 

One can also consider the supersymmetric version of this model, where there are 
also appropriate fermionic terms [11,12]. In fact, this model has also appeared as 
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a regularization of the light-cone membrane action[13]. For d — 9 it describes the 
low-energy dynamics of D-particles [14-16] and has recently been proposed as a 
matrix model description of M-theory in the large- N limit [17]. This model can 
be thought of as our model where a mass therm has been given to the Lagrange 
multiplier matrices ^e 2 A?-. Integrating out A y - (that is, solving their equation of 
motion) generates the above potential term (34) with a coefficient 1/e 2 . For e — > 0, 
corresponding to our model, the strength of the potential grows very large. There- 
fore, our model can be though of as the low-energy limit of the M-theory model, 
when the Mi are at the minimum of the potential (34), that is, they commute. In 
the super- Yang-Mills (M-theory) model there is also the Gauss law constraint 

Y^i[Mi, Mi] +^6ief = (35) 

i i 

where 6j are Majorana-Weyl spinors, the superpartners of M;. The bosonic part 
of (35) is just the "angular momentum" K. The fermionic part of (35) provides 
real representations of SU(N), consisting of (the irreducible components of) to- 
tally antisymmetric tensor products of adjoint representations (9 transforms in 
the adjoint of SU(N)). Consequently, (35) constrains the "angular momentum" 
K to be in one of these representations, depending on the fermionic sector of the 
model. Therefore we conclude that the low-energy motion of the eigenvalues of 
Mi (that is, the coordinates of the so-called DO-branes in M-theory) is the one of 
particles in the <i-dimensional spin-Calogero model (d — 9 for M-theory), where 
the spin degrees of freedom of the particles, which participate in the dynamics, 
are determined by the fermionic state of the model through (35) and (19). The 
physical implications of this fact are left for future study. 

We conclude with some remarks. Clearly, there are a lot of unanswered ques- 
tions here. Firstly, the physical meaning of the Lagrange multiplier matrices Ay 
and their role at classifying the types of solutions are yet to be understood. The 
matrix equations of motion have barely been touched in the general case, and their 
solutions are unknown. Even in the case A=constant the general solution has not 
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been fully studied. Generalizations of these models involving the Weierstrass func- 
tion potential could be sought, where instead of a matrix model one would have 
to deal with an appropriate topological model [18]. In fact, it would be interesting 
to consider what type of model would give rise to Calogero-type dynamics on a 
manifold of more general geometry and/or topology. Finally, the quantization of 
the model, being a constrained system, is a subject of investigation. 

REFERENCES 

1. F. Calogero, J. Math. Phys. 10 (1969) 2191 and 2197; 12 (1971) 419; Lett. 
Nuovo Cimento 13 (1975) 411; F. Calogero and C. Marchioro, Lett. Nuovo 
Cimento 13 (1975) 383. 

2. B. Sutherland, Phys. Rev. A4 (1971) 2019; 5 (1972) 1372; Phys. Rev. Lett. 
34 (1975) 1083. 

3. J. Moser, Adv. Math. 16 (1975) 1. 

4. D. Kazhdan, B. Kostant and S. Sternberg, Comm. Pure Appl. Math. 31 
(1978) 481. 

5. M.A. Olshanetskii and A.M. Perelomov, Phys. Rep. 71 (1981) 314; 94 (1983) 
6. 

6. J. Gibbons and T. Hermsen, Physica Dll (1984) 337. 

7. S. Wojciechowski, Phys. Lett. Alll (1985) 101. 

8. Z.N.C. Ha and F.D.M. Haldane, Phys. Rev. B46 (1992) 9359. 

9. J. Minahan and A.P. Polychronakos, Phys. Lett. B302 (1993) 265. 

10. J. Minahan and A.P. Polychronakos, Phys. Lett. B326 (1994) 288. 

11. M. Claudson and M.B. Halpern, Nucl. Phys. B250 (1985) 689. 

12. M. Baake, M. Reinicke and V. Rittenberg, J. Math. Phys. 26 (1985) 1070. 

13. B. de Witt, J. Hoppe and H. Nicolai, Nucl. Phys. B305 (1988) 545. 



11 



14. U.H. Danielsson, G. Ferretti and B. Sundborg, Int. Jour. Mod. Phys. All 
(1996) 5463. 

15. D. Kabat and P. Pouliot, Phys. Rev. Lett. 77 (1996) 1004. 

16. M.R. Douglas, D. Kabat, P. Pouliot and S.H. Shenker, Hiep-th/9608024 . 



17. T. Banks, W. Fischler, S.H. Shenker and L. Susskind, |hep-th/9610043 . 

18. A. Gorskii and N. Nekrasov, Theor. Math. Pnys. 100 (1994) 874. 



12 



